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ABSTRACT 

We measure the distribution function of rotational velocities 4){Vc) of late-type galaxies from 
the HIPASS galaxy catalogue. Previous measurements of the late-type velocity function are 
indirect, derived by converting the galaxy luminosity function using the relation between 
galaxy luminosity and rotation velocity (the Tully-Fisher relation). The advantage of HIPASS 
is that space densities and velocity widths are both derived from the same survey data. We 
find good agreement with earlier inferred measurements of 4i{Vc), but we are able to define 
the space density of objects with Vc as low as 30 kms~^. The measured velocity function 
is 'flat' (power-law slope a w —1.0) below Vc ~ lOOkms"^. We compare our results with 
predictions based on ACDM simulations and find good agreement for rotational velocities in 
excess of 100 km s~^, but at lower velocities current models over-predict the space density of 
objects. At = 30 kms~^ this discrepancy is approximately a factor 20. 
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1 INTRODUCTION 

The galaxy luminosity function is one of the best-known tools for 
comparing theoretical models of galaxy formation with the ob- 
servable universe. Modern large scale cosmological simulations, 
based on the hierarchical clustering of gravitationally interacting 
particles, need to reproduce the space density of galaxies seen in 
the observed luminosity function. However, in reality such a com- 
parison is complicated because galaxy luminosities are not a di- 
rect output of these simulations. Galaxy luminosity is non-trivially 
linked to halo mass through poorly understood processes such as 
star formation efficiency, feedback, and gas accretion. Recent semi- 
analytical models attempt to tie the A'^-body simulations to ob- 
servable quantities by applying physically motivated or empirical 
recipes for star formation, gas accretion, AGN feedback, etc. (e.g., 
Baugh et al. 2005; Bow er et al.|2006[ [PeXucia & Blaizot|2007| 
Obreschkow et al.|2009| l. These models produce a range of mea- 
surable galaxy properties, such as luminosity in different bands, gas 
masses, and star formation rates, and the distribution functions of 
these quantities can be tested against observations. Although such 
semi-analyical models are increasingly more successful in repro- 
ducing the observable universe (see review by Baugh 2006^ , the 
method inherently implies a very indirect way of comparing galaxy 
space densities. A much more direct test of the simulations can be 
made if the properties of simulated haloes are compared directly 
with observable quantities. Such a comparison is possible through 
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the galaxy circular velocity function <j>{Vc), which is defined as the 
number density of galaxies as a function of their rotational velocity. 

The galaxy circular velocity traces the total dynamical mass of 
a galaxy, which is dominated by its dark matter content. Therefore, 
the measured velocity function should match that of dark haloes in 
simulations. However, it should be realized that baryons still poten- 
tially affect the comparison in two unrelated ways. First, for surveys 
to detect the galaxies and include them in samples that are used to 
create the velocity function, the galaxies need to have a minimal 
detectable amount of baryons in the component and phase being 
searched for. Whether in an individual galaxy these baryons are 
primarily in cold gas, in old stars or young stars, may affect this 
galaxy's detectability. Second, baryon physics can significantly af- 
fect the density structure of dark matter haloes and consequently 
have an effect on its rotational velocity. Especially in the inner re- 
gions of galaxies the baryonic component can contribute signif- 
icantly to the rotational velocity and in some cases the baryons 
can dominate the dynamics in the central regions. Therefore, al- 
though the velocity function provides a much more direct tool to 
compare simulations and observations than galaxy luminosity func- 
tions, these two limitations should be kept in mind in making the 
comparison. 

There have been a number of attempts in the past to construct 
the galaxy velocity function by using published luminosity func- 
tions and then converting these into velocity number density dis- 
tributions using scaling relations such as the Tully-Fisher relation 
(for late-type galaxies) and the Faber-Jackson relation or the funda- 
mental plane (for early-type galaxies) . |Shimasaku| ( | 1 993) and |Gon-] 
|zalez et aL| l |2000^ apply this method on an array of optical surveys 
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and find consistent velocity functions (within a factor of 2), witfi a 
Schechter-like distribution, but witfi a steeper than exponential de- 
cline at the high velocities. As is stressed by Kochanek & White] 
( |2001) and demonstrated by [Sheth et al.| ( |2003| , this simple con- 
version from a luminosity function via a scaling relation can lead 
to erroneous measurements of 4'{Vc) if the scatter on the scaling 
relations is not taken into account. [Kochanek & Whitelp001[ > men- 
tion a number of other systematic problems with this method of 
measuring (t>{Vc), prime among which are the errors in morpho- 
logical classification of galaxies and calibration of the magnitude 
scale. They put special effort in reducing these systematic uncer- 
tainties and measure a non-parametric 4>{Vc). Sheth et al. j2003f 
measure 0(14) for early type galaxies from SDSS taking full ac- 
count of the scatter in the relation between luminosity and velocity, 
|Chae| j2008 1 takes a similar approach for SDSS and 2dF galax- 
ies, and Choi et arij2007^ use volume-limited SDSS subsamples to 
construct a velocity function for early types. 

In this paper we present the first measurement of the veloc- 
ity function of late-type galaxies that is not dependent on scaling 
relations. We base this measurement on the H I Parkes All Sky Sur- 
vey (HIPASS), a blind 21-cm emission line survey covering the 
whole southern sky up to declination +25°. From the region south 
of declination +2° a catalogue of 4315 Hl-selected galaxies has 
been extracted as presented in Meyer et al.| ( |2004] ) and |Zwaan et al.| 
( [2004^ . The advantage of this large galaxy sample is that the veloc- 
ity width W is one of the data products of the original survey data. 
Therefore, after doing the appropriate volume corrections and cor- 
recting the velocity widths for inclination, the velocity function can 
be measured directly without taking the step of measuring Vc sep- 
arately, or assuming a global scaling relation for all galaxies. The 
HIPASS catalogue is 99 per cent complete at a peak ux of 84 mjy 
and an integrated ux of 9.4 Jykms"^. This integrated flux limited 
translates to a minimal detectable Hi mass of 5.5 x lO^Mo at a 
distance of 5 Mpc. 

The paper is organized as follows. In section 2 we present 
the galaxy sample we use in our analysis, summarize the method 
for calculating space densities and discuss the observed velocity 
width function. Velocity functions for late-type galaxies are pre- 
sented and compared with previous inferred velocity functions in 
section 3. We also discuss the contribution from early-type galax- 
ies. In section 4 we compare our results with those of recent 
ACDM simulations. Section 5 discusses some observational bi- 
ases that could affect the evaluation of the velocity function. Fi- 
nally, we summarize our results in section 6. For distance de- 
pendent quantities we use Hq = 70 kms~^ Mpc~^ and define 
hro = Ho/{7Q kms"' Mpc^^) = 1. 



2 THE HIPASS VELOCITY WIDTH FUNCTION 
2.1 Calculating space densities 
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In jZwaan et al.| ( (2003| l we derived a bivariate stepwise maximum 
likelihood (2DSWML) technique, which solves for the space den- 
sity of objects as a function of H I mass and velocity width simul- 
taneously. We specifically designed this method for a peak-flux- 
limited sample, but in [Zwaan et al.| ( [2004[ > we demonstrated how 
the method can be adjusted to produce reliable space densities for 
samples constructed using different selection criteria. In that paper 
we also showed how to use the 2DSWML technique to calculate 
for each galaxy individually an effective volume Vcs. These values 
are now the maximum likelihood equivalents of Knax for the stan- 
dard El/Knax method i jSchmidt 1968 1 and summing the values of 
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Figure 1. The space density of local gas-rich galaxies as a function of 
their observed profile widths corrected for relativistic broadening and in- 
strumental broadening. The short-dashed line is a modified Schechter fit to 
the data points (see text). The solid line is the distribution function based 
on the |Obreschkow et al.H2009) simulations, while the long-dashed line 
represents this sample limited to late-type galaxies only. 



1 /VcS in bins of H I mass gives the original 2DSWML solution for 
the H I mass function. The normalization of space densities is re- 
tained by equating the mean galaxy density n to the integral over 
the Hi mass function. After this normalization is applied to T4ff, 
these values represent true volumes. Any distribution function of 
a parameter measured for the HIPASS sample can now be deter- 
mined by a binned summation of l/VcS- We refer to Zwaan et al 
(2004 J for a more detailed description of the method. 



2.2 The velocity width function 

We use the velocity widths W measured at 50 per cent of the peak 
flux density of the H I profile. In principle, for high signal to noise 
data, the widths measured at 20 per cent of the peak flux density 
will provide more direct measurements of the full rotation, but we 
prefer the 50 per cent values because these latter values are less 
sensitive to noise in the 21-cm spectra. [Meyer et li] ( |2008] l have 
shown that the scatter in the TuUy-Fisher relation is smaller when 
using 50 per cent measurements instead of those measured at 20 
per cent. The observed line widths are corrected to their rest-frame 
values (relativistic effect), instrumental broadening and turbulent 
motion, all according to standard procedures. A full description of 
the corrections applied to the observed velocity widths is described 
in |Meyer et al.| p008). Henceforth, we refer to the velocity width 
corrected for instrumental and relativistic broadening as W. Note 
that for the calculation of rotational velocities an additional correc- 
tion for turbulent broadening is applied. 

Fig[T] plots the distribution function of H I profiles widths, 
(piW). This function can be readily calculated by summing 1/Vca 
in bins of W. Error bars represent uncertainties in the measure- 
ments determined by quadratically summing 1/VcS within each 
bin. This estimate gives an indication of the pure statistical uncer- 
tainty, which mostly ignores effects due to local large scale struc- 
ture. All 4315 HIPASS detections presented in Meyer et al.H2004| l 
are represented in this diagram as no optical identifications are re- 
quired for constructing this function. The error on W m easure- 
ments is aw = 7.5kms~^, independent of profile width i Zwaan 
let al.|2004| >. 
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To provide a parametric description of this distribution, we fit 
it with a modified Schechter function of the form 



W . 



<t){W)dW = MivrT exp[-(-^) 



T(q//3) W 



(1) 



where F is the gamma-function. In the next section we give a phys- 
ical motivation for this parameterization of a velocity function. We 
find a good fit to the measured data points with the following val- 
ues: = 0.67/i^oMpc~^, logH^./kms"^ = 2.21, a = 0.10, 
and/3 = 1.56. 

Although no true physical parameter is represented in this di- 
agram, it does present a straightforward comparison to simulated 
galaxy samples. Since it does not depend on cross-matching with 
optical catalogues, it shows the distribution function of observed 
profile widths of a purely H I-selected sample. This function is rel- 
evant for calculating the detection efficiency of future Hi 21 -cm 
surveys, as the ability to separate 21 -cm signals from the noise is a 
strong function of profile width. 

An interesting comparison is with the work of |Obreschkow| 
|et al.| ( [2009} , who calculate cold gas properties for the semi-analytic 
galaxy catalogues of |De Lucia & Blaizot| ( 2007^ based on the A'^- 
body simulations of cold dark matter by Springel et al.| ( |2005| l. 
Velocity widths are assigned to simulated galaxies by making as- 
sumptions on the rotation curves and gas distributions. We refer to 
|Obreschkow et aL| ( [2009[ > for details on how velocity widths are 
calculated. The solid line in Fig. [T] shows the distribution function 
of W from these simulations. Below W — 100 kms^^ we draw 
a dashed line to indicate that the simulations are incomplete in that 
range due to the limited mass resolution. For the W measurements 
from |Obreschkow et al.| {2009} we used the raw W50 values multi- 
plied by sini. In the range 100 kms~^ < W < 400 kms~^ the 
HIPASS points match the |Obreschkow et al.| ('2009l points fairly 
well, although the simulations over-predict the space density of 
galaxies just above the mass resolution limit. Toward higher ve- 
locities the observations show an exponential decline, whereas the 
observations display a power-law distributions. Clearly, the simu- 
lations predict too many very broad profile (W > 400 kms~^) 
galaxies compared to what HIPASS detected. The H I mass func- 
tion of the |Obreschkow et al.|p009i catalogue also shows an over- 
production of very massive galaxies (their Fig. 3), but the differ- 
ence in our Fig[T]is much larger. To investigate this discrepancy 
further, we also plot in the same figure the simulated distribution of 
late-type galaxies only, and find that the match the observed distri- 
bution (without any morphological type cut applied) is much better. 
Therefore, it appears that the |Obreschkow et al.| ( |2009l l simulations 
overproduce a population of gas-bearing massive early-type galax- 
ies that is not observed in HIPASS. 



3 A TYPE-DEPENDENT ROTATIONAL VELOCITY 
FUNCTION 

3.1 A velocity function for late-type galaxies 

The next step is to determine the distribution of true rotational 
velocities. These 14 values we calculated by applying an incli- 
nation correction to the projected velocity spread of the Hi gas: 
Vc = W{2sm i)"^ Henceforth, we refer to the rotational veloc- 
ity Vc as the velocity width corrected for turbulence, instrumental 
broadening, inclination and corrected to rest-frame profile widths, 
as detailed in |Meyer et al.|j2008) . We note that K should be re- 
garded as the velocity corresponding to the maximum of the ro- 
tation curve as this is essentially the velocity that is probed by 



integrated Hi velocity profiles. In principle, more accurate mea- 
surements of galaxy halo masses could be achieved my making use 
of Vflat, the velocity corresponding to the flat part of the rotation 
curve, but high resolution 21cm synthesis observations are required 
to make such measurements . [\^rheij en| pOO 1 | l discusses the use of 
Vc and Vflat in the TuUy-Fisher relation and their relation to halo 
mass. 

The required inclinations cannot be measured from the 
HIPASS data directly because the Parkes Telescope beam at 21- 



cm of 13' (FWHM) does not spatially resolve most galaxies. Doyle 



|et al.| ( (2005] ) searched for optical counterparts to the HIPASS detec- 
tions and found unique optical identifications for 2646 (62 per cent) 
galaxies. For these galaxies [Doyle et aL] j2005] l list 6j-band axis- 
ratios measured from the SuperCOSMOS data. Inclinations are cal- 
culated from the axis ratios q and corrected for internal thickness 
of the disks go using the relation 



1 



(2) 



where we adopt go — 0.2. It has been argued that go should de- 
pend on galaxy morphology (Giovanelli et al. 1997 ) but we fol- 
low the recommendation of ,Tully et al.,(2009) that a consistent use 
of go is important for the avoidance of systematic errors. We note 
that [Meyer et al.H2008^ found that inclinations derived from the 
bj measurements resulted in lower TuUy-Fisher scatter that those 
from the 2MASS J-band measurement. 

Selecting galaxies on the basis of their H I content naturally re- 
sults in a sample weighted toward gas-rich galaxies. However, this 
does not imply that early-type galaxies (SO/E) are missing from the 
catalogue. In fact, some recent studies have shown that 70 per cent 
of early-type galaxies have detectable amounts of H I with masses 
exceeding a few times 10*' Mq (e.g., Morganti et al.|2006 



Gost- 



[erloo et al.|2007[ l. The H I morphologies range from small clouds 
to large, regular low column density disks. Furthermore, in |Zwaan] 
[etaL[f2003 l we have measured the H I mass function of early-type 
galaxies, based on an initial sample of 1000 bright HIPASS galax- 
ies, indicating that early-types are present in 21 -cm survey cata- 
logues. However, in the analysis of the HIPASS data we wish to 
exclude the early-type galaxies for two reasons. First, even though 
a large fraction of early-types has detectable H I, we are certainly 
missing a significant fraction of those early-types that are gas free 
or sufficiently gas poor to be absent from the HIPASS sample. Sec- 
ondly, Hi measurements do not provide reliable velocity disper- 
sions for early-type galaxies. As already stated, early-type H I mor- 
phologies are not predominantly ordered disks that trace the full po- 
tential of the galaxies. This fact is also supported by the observation 
that H I masses do not correlate with early-type luminosity or mass 
l |Serra et aT. 2008 1. Optical surveys such as the SDSS, completed 
with optical spectroscopy are much better equipped to measure the 
distribution function of early-type velocity dispersions ( [Sheth et aT] 
[2003[|Choi et al.|2007[ [Chae|2008| l. To cull the early-type galax- 
ies from the sample we rely on the morphological classifications 
as listed by [Doyle et al.| l [2005^ . This takes out 1 1 per cent of the 
galaxy sample. 

We also exclude all galaxies with measured inclinations i < 
45° to avoid introducing large uncertain inclination corrections that 
may affect the shape of the velocity function. This amounts to 30 
per cent of the sample. To account for the excluded galaxies, the 
final velocity function is corrected with a factor 1/(1 — /) where / 
is the fraction of excluded galaxies. 

In Fig |2] the measured velocity function is shown as solid 
circles, whereas the open symbols indicate the effect of chang- 
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Figure 2. Top panel: The HIPASS late-type velocity function 4i(Vc). The 
solid circles show (/'(Vc) for morphological types Sa and later and inclina- 
tions i > 45°. The other symbols show the effects of choosing different 
inclination cut-offs and including early-type galaxies. The solid line rep- 
resents the best-fit modified Schechter function (see text). Bottom panel: 
Histogram of rotational velocities for our sample of HIPASS galaxies with 
i > 45° and types Sa and later 



ing the inclination cut-off from 45 to 35 and 55 degrees. The ef- 
fect of changing the inclination cut-off is not noticeable for K 
above approximately 100 kms~^, but at lower Vc the changing 
cut introduces some scatter. There is a weak trend of space den- 
sities dropping with a higher inclination cut, and increasing with 
a lower inclination cut, although the error bars are too large to 
make a conclusive statement: the mean space density of objects 
with Vc < 80kms~^ changes less than la when the inclination 
cut is varied from 45 to 35 or 55 degrees. 

Also in Fig|2] we show the effect of ignoring the morpho- 
logical type selection by plotting the velocity function for all 
morphological types. The effect is only detectable beyond K > 
200 kms"\but even there early-type galaxies make a very minor 
contribution to the space densities. 

Similar to the distribution function of W , we fit the rotational 
velocity function with a modified Schechter function of the form 



0(Vc)dyc = 0,( 



(3) 



where is the space density of galaxies near the 'knee' of the 
distribution, K, is the characteristic rotational velocity, a is the 
low-velocity power-law index, and /3 is the high-velocity power- 
law index. The latter parameter l3 enters the equation when the tra- 
ditional Schechter luminosity function is transformed to a veloc- 
ity function using a relation between luminosity L and rotational 
velocity, e.g., the TuUy-Fisher relation: L oc Vc'^ . In reality, the 
value of P in the parametrized velocity function is not equal to the 
Tully-Fisher power law because the scatter in the Tully-Fisher re- 
lation affects the way in which the luminosity function is trans- 
formed to <j){Vc), and because the Tully-Fisher power law changes 
as a function of luminosity and morphological type. These effects 
have been described in detail by e.g., |Sheth et al.| ( |2003^ . A tra- 
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Figure 3. The HIPASS late-type velocity function 4){Vc) as indicated by 
the circles with error bars. The solid curve is a modified Schechter function 
fit. The two dashed curves show the late-type velocity function based on the 
2dF data (upper curve) and SDSS data (lower curve) as calculated by |Chae| 
(2008) . 



ditional Schechter function would provide a very poor fit to the 
measured velocity function because the high end exponential drop 
off does not decline fast enough to fit the data points. We find a 
satisfactory fit to the data with the following parameters: 0* — 
(0.061 ± 0.027) /^?oMpc-^ log K./kms'^ = 2.06 ± 0.09, 
a = 0.66 ± 0.47, and P = 2.10 ± 0.37. It should be noted that 
errors on individual parameters are very large because of the strong 
covariance between all fitting parameters. In particular, the dimen- 
sionless covariance between logK and P is r(logK,/3) ~ 0.99 
and that between 0* and a is r{<l), a) — —0.97. The lowest covari- 
ance value is r((;/>, , P) — 0.78. By just fitting the low velocity slope 
(K < 100kms"^)of the distribution, we find a = -0.90±0.25. 

Fig.|3]reproduces the HIPASS velocity function and compares 
our result with recent estimates of the late-type velocity function 
based on the 2dF and SDSS galaxy redshift surveys from |Chae| 
(2008 ). In that study, published luminosity functions were adopted 
for the SDSS late-type galaxies from jChoi et al.|p007^ and for 2dF 
from |Croton et al.|(2005^ . Next, they Monte Carlo simulate a large 
sample of galaxies, based on these luminosity functions and the 
Tully-Fisher relation, taking into account the intrinsic scatter. For 
the Tully-Fisher relation, the measurement by Pizagno et al. ( 20()7jl 
was adopted, which is based on Ha rotation curves of a sample of 
162 SDSS galaxies. The velocity functions are found by summing 
the simulated galaxies in bins of Vc- In Fig. [5] the upper dashed 
curve shows the result based on the 2dF luminosity function and 
the lower dashed curve is based on SDSS. Discrepancies between 
the two function are indicative of the systematic errors introduced 
by converting galaxy luminosities between different wave bands 
that are used in the luminosity function and Tully-Fisher relation 
measurements, different definition of morphological type (based 
on spectroscopy or colour), and differences in the 2dF and SDSS 
galaxy luminosity functions (as discussed in e.g., |Liske et al.|2003] ). 

Given the very different galaxy samples and methods for cal- 
culating the velocity function, we find reasonable agreement be- 
tween the directly measured HIPASS velocity function and the in- 
ferred 2dF and SDSS velocity functions. At lower velocities, we 
find slightly higher space densities compared to the indirect mea- 
surements and at higher velocities we find slightly lower values. 
The most significant difference between our measurement and pre- 
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vious estimates is that we can determine ttie space densities of 
galaxies down to much lower values of T4: the HIPASS veloc- 
ity function goes down to 14 « 30kms~^, whereas the previous 
measurements only reach K ~ 70kms~^. The measured velocity 
function is remarkably 'flat' (power-law slope ~ — 1) down to the 
lowest rotational velocities. 



3.2 The contribution from early-type galaxies 

Because 2 1 -cm emission is not a good tracer of the gravitational po- 
tential of early-type galaxies, we have so far excluded these galax- 
ies from our analysis. It would, however, be interesting to derive 
the total local universe velocity function of all galaxy types. To 
put early and late-type galaxies on the same K scale, we adopt 
a conversion from velocity dispersion a to rotational velocity: 
K = ^/2a. This conversion is based on the assumption that the 
mass distribution of an elliptical galaxy can de modeled as a singu- 
lar isothermal sphere. |Desai et"aLl ( |2004^ test this assumption us- 
ing literature measurements of velocity dispersion and optical cir- 
cular velocities in the flat regions of the rotation curves and find 
Vc ~ 1.54(7, illustrating the uncertainty in estimating K values for 
early-types. 

The most recent estimates of early-type velocity functions are 
those by |Chae| POOS'!, based on 2dF and SDSS using a method 
similar to that described for the late-type galaxies. Fig. |4] shows 
these measurements as indicated by a shaded thick curve, of which 
the upper boundary corresponds to the calculation based on 2dF 
data and the lower boundary corresponds to the SDSS. Other re- 
cent calculations of the SDSS velocity function have been pre- 
sented bylSheth et al.|^2003|l and |Choi et"ar] ( |2007| l, the latter be- 
ing based on direct measurements of velocity dispersions. How- 
ever, this study may suffer from incompleteness at low velocity 
dispersions. As shown by |ChaeH2008| >, the different space density 
measurements of early-type galaxies around the peak of the distri- 
bution (K ~ SOOkms"^) a gree very well, but below that point 
the slopes of </'(14) vary between the different studies. Issues with 
completeness, both related to incomplete dispersion measurements 
and missing low luminosity early-types, are probably causing the 
discrepancies. See also [Mitchell et ak] ( |2005^ for a discussion on 
inferred and measured early-type velocity functions. 

The solid curve in Fig.|4]shows the summed velocity function 
for all morphological types, adding the HIPASS late-type curve and 
the 2dF/SDSS-based early type curve of Chae ( 2008 1. For compari- 
son, we also show the velocity function from Gonzalez et al. ( 2 000) 
based on SSRS2 data. Their result is in good agreement with later 
calculations by |Kochanek & White| ( [200T] l, although these authors 
do not give a parametrized description of their velocity function. 
The overall consistency between the HIPASS plus 2dF/SDSS ve- 
locity function with the earlier work of [Gonzalez et al.| (2000) is 
reassuring. Again, HIPASS adds important data points at the low 
Vc end of the velocity function. 

Fig.|4]illustrates that the galaxy velocity function is dominated 
by late-type galaxies below K ~ 200 km s^^, whereas early-type 
galaxies take over above that point. It should be noted that the early- 
type measurements are based on samples having relatively high 
luminosity cut-offs. Dwarf elliptical (dE) galaxies are therefore 
largely missing from the analysis. In field regions late type-spirals 
and irregular (Irr) galaxies dominate the space densities at the low 
mass end. This is both apparent from optical surveys ( Bin ggeli et al.| 
[19881 and H I surveys ( Zwaan et al.|2003l l. Dwarf ellipticals only 
become significant contributors to the space density in cluster envi- 
ronments. The evidence from H I mass function studies is that 21cm 
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Figure 4. The points with error bars show the HIPASS late-type velocity 
function with a modified Schechter function fit indicated by a short-dashed 
line. The shaded curve represents the early-type velocity function derived 
from 2dF data (upper boundary) and SDSS data (lower boundary) as esti- 
mated by [ChaeH2008) . The solid curve is the summed velocity function for 
all morphological types. The [Gonzalez et aT]j2000) all-types result is rep- 
resented by a long-dashed line. Late-type galaxies dominate the statistics 
below Vc Si 200kms"i. 



surveys fully recover the the faint-end slope of the luminosity func- 
tion ( [Zwaan et al.|200l| l. Therefore the contribution from gas-less 
early-type dwarfs is marginal in the volume sampled by HIPASS. 
[pesai et al. ( 2004 ) construct velocity functions for cluster galaxies, 
and find that their function is steeper than that in the field, which 
is also attributed to a contribution of low mass early-type galaxies, 
which are relatively scarce in a large volume of space. 



4 COMPARISON WITH THEORY 

The cold dark matter model of galaxy formation with a cosmolog- 
ical constant (ACDM) makes robust predictions for the mass func- 
tion of dark matter haloes. Cosmological simulations determine the 
mass function slope with high accuracy, giving ip{M) oc M~^'^ 
(e.g.,[Sheth e t al. 2001 , Sheth & Tormen 2002, Gi ocoli et al.pOOSl 
[Tinker et a l 2008; Boyla n-Kolchin et al.| 2009 1. To convert this 
function into the distribution function of individual galaxies, halo 
occupation distribution (HOD) has to be applied ( [Peacock & Smith| 
[2000[[Cooray & Sheth|2002t . The number density of galaxies may 
be higher than that of dark matter haloes since each dark matter 
halo can contain a number of galaxies. 

ACDM models generally predict that the total dynamical 
galaxy mass Mdyn relates to rotational velocity K as Afdyn oc Vc*^. 
In our definition of rotational velocity, K is comparable to the max- 
imum rotational velocity, for which definition [Bullock et al.[(200T] l 
found that P ~ 3.4. Applying this conversion to the power-law 
mass function slope of qm = —1.9, we can estimate the power- 
law slope a of the velocity function 



<t){Vc)dVc = i>{M)dM 



(4) 
(5) 
(6) 



finding that power-law slope of <j!>(K) should be approximately 
a — —4.0. Of course, this simple estimate ignores many subtleties 
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of HOD and baryon physics, but is does illustrate that theory prob- 
ably predicts a low-K slope much steeper than q ~ — 1, which we 
find for the HIPASS velocity function. 

The clearest manifestation of this problem with ACDM is 
the so-called substructure problem: around the Milky Way Galaxy 
ACDM predict hundreds of low-mass satellites, although only an 
order of magnitude fewer objects are detected l |Klypin et al.|1999l 
[Moore et aH1999| . Recent surveys identify a growing number of 
satellites around the Galaxy and M31 (see jZucker et al.|2006[|Ibata| 
let al.|2007| for some examples), but still not in sufficient quanti- 
ties to alleviate the problem. By accounting for baryonic physics, 
|Simon & Geha| l (2007[ > manage to reduce the missing fraction of 
satellites to approximately a factor of 3, but baryon physics are 
extremely complicated and uncertain in these dense environments 
to the extent that reconciliation between theory and observations 
is presently too model dependent. We note however that |Koposov| 
|et al.| l |2009^ demonstrate that models with a constant stellar mass 
fraction and strong suppression of star formation before reioniza- 
tion in haloes with K < 10km can produce the observed satel- 
lite counts. 

The range of masses that we are probing with HIPASS is 
clearly higher than that of the Milky Way satellites. H I-selected 
galaxy samples are inherently composed of fairly isolated systems 
as 21-cm surveys are not very efficient in probing high density re- 
gions ( jWaugh et al.|2002l l. The galaxies that we probe in the tail 
of the velocity function are predominantly dwarf irregular galaxies 
and low mass spirals (cf. |Zwaan et al.^2003j . 

To make a direct comparison with simulations, we need a 14 
distribution function of dark matter haloes. [Blanton et al.| j2008j 
compare an inferred velocity function of isolated dwarf galaxies 
with isolated haloes in the TV-body, pure dark matter simulations 
of |Kravtsov et al.| ( [2004[ >. Due to the first order similarity of the 
observed data sets, we also use this Kravtsov based sample as our 
dark matter halo reference distribution. 

Unfortunately, Blanton et al.|p008] l present cumulative dis- 
tribution functions of K, with makes quantitative comparisons 
very difficult. We differentiate their cumulative distribution func- 
tion of dark matter haloes to derive (j)(Vc) and show the result 
in Fig |5] as a dotted line. Interestingly, the space densities calcu- 
lated from simulations and HIPASS match very well in the range 
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< Vc < 200 km s . To some level, this is a co- 



incidence because the definition of isolation for the dark matter 
haloes is not matched to that of the HIPASS catalogue. |Blanton| 
|et al . ( 2008 ) show that the space densities in this K range increase 
somewhat if a more relaxed definition of isolation is applied. Of 
more significance is the marked difference in the slope of the veloc- 
ity function below Vc ~ lOOkms"^, where the simulations pre- 
dict a slope of a ~ —3, compared to our observation of a ~ —1. 
Simulations clearly predict many more haloes than are detected by 
HIPASS. 

There is one important caveat in this comparison, which is 
the mass resolution of the cosmological simulation. The ( [Kravtsovl 
|et al. 2004 1 simulations used by [Blanton et"ar] ([2008) to calcu- 
late velocity functions, reliably determine space densities of haloes 
down to K ~ 70kms~^. Below that value, we have used the 
power-law extrapolation of |Blanton et al.| ( [2008) . 

Another interesting comparison with our observations is that 
with the recent results of .Obreschkow et al.i (20091. These au- 



thors have reprocessed the semi-analytic simulations of ]De Lucia &| 
|Blaizot| ( [2007[ l, which, in turn, are based on the 'milleimium run', a 
cold dark matter (CDM) A^-body si mulation of a (714 Mpc/h 7o)^ 
volume, containing 3 x 10^ galaxies. Obreschkow et al. ( 2009 1 use 
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Figure 5. The HIPASS velocity function (filled circles with solid line) com- 
pared with the velocity function from a selection of simulations. The blue 
dashed line shows the simulations of |Obresclikow et ar]j2009f for late-type 
galaxies. The green short dash-dotted and long dash-dotted lines are the 
CDM and WDM simulations from |Zavala et al.H2009^ for all galaxy types. 
For comparison the observed velocity function for all galaxy types from 
Fig|4]is reproduced as the black dashed line. The red dotted line shows the 
dark matter haloes from ,Blanton et ahljlOOS) . 



a set of physical prescriptions to assign H I and molecular hydro- 
gen (H2) to these galaxies and are able to reproduce many of the 
observed gas properties of local galaxies. The simulated H I mass 
function is well matched to the one measured from HIPASS for H I 
masses in excess of A/hi ~ 2 x IO^Mq. The scatter in the ob- 
served and simulated A/hi — W50 relation is very large (their Fig. 
13), but this A/hi limit approximately translates into a rotational 
velocity limit of K ~ lOOkm s"'^, implying that space densities of 
galaxies below K ~ 100km are not reliably predicted. 

To estimate K values from the |Obreschkow et al.| ( |2009l l cata- 
logue, we use their measurements of ui^'^j , which is the line width 
between the two peaks of the double homed Hi profile. We use 



peak 



+ 12)/2 to derive rotational velocities, where the 
correction of 12 km is applied to correct for profile smoothing 
due to the velocity dispersion of the gas. That is, w^'^^j^ slightly un- 
derestimates the true K as can be seen in Fig. 1 1 of |Obreschkow| 
[etaL] ( |2009) . 

The thick dashed line in Fig. |5] shows the resulting velocity 
function, plotting the result for late-type galaxies to enable a direct 
comparison with the HIPASS results. The agreement in the shape 
of ^(Vc) above K ~ 100 km s^^ is striking, only above K — 
300kms~^ do the simulations over-predict the observed number 
density. Unfortunately, at the low velocity end we cannot make a 
detailed comparison because the mass resolution of the millennium 
run is not sufficiently high to predict space densities in this regime. 

Higher resolution simulations are available from the recent 
work of [Zavala et al.| ( [2009 l. These authors use constrained cold 
and warm dark matter (WDM) simulations of a {29h^Q Mpc)"^ 
volume centered around the Local Group. The WDM simulations 
are based on a power spectrum using a fitting function that ap- 
proximates the transfer function for a thermal WDM particle with 
mwDM ~ IkeV. Their velocity functions are reproduced by short 
dash-dotted and long dash-dotted lines for CDM and WDM, re- 
spectively. Although their simulations assume that all galaxies are 
disk galaxies, they do not differentiate specifically between early or 
late type galaxies. We reproduce the velocity function for all galaxy 
types from Fig|4]as a long dashed line to compare the [Zavala et al.| 
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([2OO9J results with. The simul ated curves are puq^ orted to be com- 
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(2009 



also test the ef- 



fect of supernova feedback on the shape of the velocity function 
and conclude that it only marginally reduces the space densities at 
the low Vc end. 

For a comparison with their simulations, [Zavala et al.H2009^ 
use preliminary results from the ALFALFA survey l |Giovanelli| 
|et al . 2005 ), concentrating on publicly released data bases in the 
Virgo direction region l^Giovanelli et al.|2007^ and the anti- Virgo 
direction region jSaintonge et al.|2008^ . In total, the comparison is 
based on a small sample of 201 galaxies, 186 of which are in a high 
density region and the remaining 15 in a low density region. Only 
galaxies out to a distance of 29h^Q Mpc are considered. Since the 
preliminary ALFALFA results are based on a much smaller galaxy 
sample, probing a smaller volume than HIPASS and concentrate 
only on a clear overdensity and an underdensity, we do not attempt 
to make a direct comparison with our results. We note, however, 
that the shape of the ALFALFA 0(K) as presented by |Zavala et al 



l |2009| l is qualitatively in good agreement with our results. 

Comparing our HIPASS measurements with the [Zavala et al.| 
( [2009) simulations, we essentially reach similar conclusions: above 
Vc = 100 km s'^ the CDM predictions are in good agreement with 
the observations. Below K = lOOkms"^ CDM clearly predicts 
too many low mass haloes. Similar to the [Blanton et al.[ ^2008) 
results, CDM predicts a steep low velocity slope as estimated by 
Eq.|6] The WDM simulations predict a fairly flat velocity function, 
in much better agreement with the observations. 



5 BIASES 

In making comparisons between our HIPASS results and simula- 
tions, there are a number of biases that should be discussed here. 
First, an important issue is whether the 21 -cm velocity widths truly 
probe the full range of the galaxy velocity curve. If the detectable 
H I is confined to the central regions of the galaxies, the H I rotation 
curve might not measure the full rotational velocity of the dark mat- 
ter halo (e.g., [Swaters et al.|2009[>. Inspecting the dwarf galaxy rota- 
tion curves of Swaters et al. (2009 ), it can be seen that most galaxies 
do just reach a flat part, but there are notable exceptions, especially 
those with rotational velocities lower than Vc « 60kms~^. This 
implies that a fraction of our lower Vc measurements probably un- 
derestimates the true rotational velocity of the dark matter haloes. 
As an extreme test, we randomly multiply half our K measure- 
ments below 60kms~^ with a factor 1.5 and construct (l>{Vc) us- 
ing this new estimates. Fig|6]shows the result, where filled symbols 
show the original space densities and open circles represent those 
after making the correction. Not surprisingly, the space density of 
objects with Vc just below 90 kms~^ is boosted after making the 
correction. In effect, the discrepancy between ACDM theory and 
our observations is slightly reduced just below the knee of the dis- 
tribution, but not significantly at lower K- The effect that the space 
density drops below 50 kms^^ is a result of the fact that many Vc 
below that value shift to higher Vc bins, but no very low velocity 
systems replace them because they do not exist in our sample. 

A second concern is whether gas velocity dispersions in low 
mass galaxies should be taken into account in calculating rotational 
velocities that measure the full potential. In very low mass galaxies, 
the gas velocity dispersions can be comparable in magnitude to the 
measured rotational velocities, adding a significant pressure sup- 
port term compared to rotational support. If H I synthesis maps are 
available, asymmetric drift correction can be made that correct for 
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Figure 6. Testing the effect of possible biases on the shape of the rotational 
velocity function. The solid circles and line show the velocity function as 
presented in Fig. |2] The open circles demonstrate the effect of increasing 
the rotational velocity of low mass galaxies for which the flat part of the ro- 
tation curve might not be reached. The open triangles show how the velocity 
function changes if velocity widths are corrected for pressure support. 



the pressure support. |Begum et al.[ ( [2008 a'bl) study a sample of very 
low mass galaxies with rotational velocities as low as 10 kms~^. 
For each galaxy individually, they correct for pressure support and 
find significant corrections to Vc- We find that their corrections can 
be fitted as AK = 130/K, which we can apply to our data to make 
a rough correction for pressure support. The triangles in Fig|5]show 
combined effect of this correction and the correction described in 
the previous paragraph. We find that the pressure correction has no 
measurable effect on (j>{Vc). 

Finally, we discuss how much cosmic variance can affect our 
measurements. In Fig. |7]we display the velocity function for four 
different quadrants of the southern sky. For each quadrant, we sim- 
ply use the El/Vcff method as described in section 2, where the 
Vcfi values are the maximum likelihood equivalents of the effec- 
tive search volumes for each galaxy, calculated from the com- 
plete HIPASS sample. Also shown are the least-squares modified 
Schechter fits to each set of data points. Although the variation be- 
tween the data points for the different quadrants is approximately 
a factor 2 at the high K end and even larger at the low Vc tail, the 
fitted functions show good consistency between the quadrants. The 
only deviant velocity function is that of the fourth quadrant, which 
is clearly an underdense region. Despite a factor two lower over- 
all normalization, the shape of the velocity function of the fourth 
quadrant is not significantly different from those of the other three 
quadrants. We conclude that our main result that the low Vc slope of 
the velocity function is 'flat' is unchanged if we subdivide our sam- 
ple into smaller regions on the sky. Cosmic variance in the HIPASS 
survey volume is discussed in more detail in [Zwaan et aL[p003) . 



6 SUMMARY AND DISCUSSION 

We have used the HIPASS sample of H I-selected galaxies to cal- 
culate the circular velocity function (t>(Vc) of late-type galaxies. 
Optical identifications of the HIPASS galaxies were used to deter- 
mine inclinations and morphological types. Our determination of 
the rotational velocity function is the first that uses a direct method, 
where the same galaxy sample is used for measurements of Vc and 
the evaluation of space densities. Previous estimates of the (j>{Vc) 
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Figure 7. The effect of cosmic variance on the circular velocity function. 
The four different symbols show the velocity function for four different 
quadrants of the southern sky, and dotted lines show modified Schechter 
fits. The solid line represents the fit to the full HIPASS sample. 



were dependent on a conversion of the galaxy luminosity func- 
tion and a relation between luminosity and rotational velocity (the 
Tully-Fisher relation). In these analyses different data sets were 
used for the calculation of the luminosity function and the Tully- 
Fisher relation. 

We find good agreement between our velocity function and 
inferred velocity function from the literature, but we are able to 
measure space densities down to Vc — 30kms^^, a factor of two 
lower than previous estimates. A standard Schechter function is not 
a good fit to 4i{Vc) because space densities decline too rapidly at 
high velocities compared to an exponential drop-off. Instead, we 
use a modified Schechter function, which provides an excellent fit 
to the data points. In order to find a velocity function for all galaxy 
types, we add the early-type measurements from |Chae|j2008^ . 

As stated in the introduction, the velocity function is one of 
the most direct tests of galaxy formation models. Our velocity func- 
tion is compared with those derived from ACDM galaxy formation 
simulations and we find a good agreement in the space density of 
objects with circular velocities > 100 kms~^. Below this value 
ACDM predicts too many objects, reaching a discrepancy of a fac- 
tor 20 at Vc — 30kms^^. WDM simulations are constructed to 
reduce the power on small scales and hence produce velocity func- 
tions with a flatter low-end slope, which provide a better fit to our 
data points. [Primack|p009 ^ argues that observational limits on the 
mass of the warm dark matter particle are already very strict. The 
circular velocity function essentially tests cosmological models on 
scales of approximately 1 Mpc ( [Blanton et al.|2008l l. Similar scales 
are probed by the Lya forest power spectrum and it is found that 
the mass of the WDM particle should be more than 4keV I Viel et af] 
2008^ (or more than 1.7keV according to a more recent analysis by 



Boyarsky et"ar|p009|>) 



Blanton et al. (2008} use a sample of 12 low mass isolated 



dwarf galaxies and an independent measurement of the luminosity 
function to compare with ACDM predictions and claim cons istency 



Zavala 



down to Vc ~ 50krns^^. Our analysis, and also that of 
|et al.| ( [2009) , indicates that the inconsistency between observations 
and simulations is significant over a larger range in velocities, up 
to K ~ 100 kms"'^. We argue that our constraints are firmer than 
those of |Blanton et al.|p008[ l because our HIPASS measurement is 



based on ~ 50 times more galaxies below K — 100 km and is 
not reliant on optical luminosity functions. 

Apart from adopting WDM, which other factors could be re- 
sponsible for the mismatch between observed and simulated ro- 
tational velocity functions? In section 5 we tested whether obser- 
vational biases, such as dwarf galaxy rotation curves not probing 
the full potential or ignoring pressure terms could have an effect 
on (t>(Vc) and concluded that the effect was insignificant. There 
are however astrophysical effects that could conceivably affect the 
shape of the measured <^(14). 

Our analysis is based on the assumption that low mass galaxies 
contain a detectable amount of neutral gas. Any gas-free low mass 
field galaxies would escape inclusion in the HIPASS sample. |Geha| 
[et al.,^2006^ obtained Green Bank Telescope 21cm observations of 
101 low mass field galaxies and detected H I emission in 88 per cent 
of those. It therefore seems unlikely that many low mass systems 
are missing from our analysis. A correction to the low K end of a 
factor 1/0.88 would not change the slope significantly. 

UV background photoionization of primordial gas could im- 
pact the formation of low mass objects (e.g., |Efstathiou|1992[|BabuT| 
|& Rees|I99 2). Recent calculations do not reach consistent conclu- 
sions on the magnitude of this effect and the characteristic K below 
which the effect becomes important. For example, variations in the 
reionization history, the inclusion of self-shielding in the calcula- 
tions, and the link with SN feedback all have important effects (e.g., 
Susa & Umemura 2004 , Hoeft et al. 2006 ; Mesinger & Dijkstra| 
2008, Sawala et al. 2009 , Pawlik & Schaye 2009 , and references 



therein). At present it remains unclear whether the suppression of 
low mass system formation by the UV background is sufficient to 
explain the observed local velocity function. Similarly, the effect 
of SN feedback might be important. SN feedback implies chemi- 
cal enrichment, heating and entraining of surrounding gas, reduc- 
ing the gas fractions of low mass systems (e.g. |Dekel & Silk|I986[ 
[White & Frenk|I99T] l. The treatment of SN feedback in simulations 
of galaxy formation is very complicated and, like UV photoion- 
ization, the importance of the effect is still a topic of discussion 
([Scannap ieco et ari|2006| l. It should be noted though that |Zav^a| 
|et al. ( 2009| conclude that the SN feedback has only a minor effect 
on the z = circular velocity function. 

Future blind Hi surveys with instruments such as the Aus- 
tralian SKA Pathfinder (ASKAP) will probe much larger volumes 
than HIPASS with a lower noise level and much higher spatial res- 
olution. For instance, a one year survey covering 75 per cent of the 
sky will find approximately 5 x 10^ galaxies, 3 x 10* of which 
will be spatially resolved, thus allowing the measurement of rota- 
tion curves and hence of more accurate rotational velocities. This, 
in combination with the much larger number of galaxies, a higher 
sensitivity to low mass galaxies, and the ability to study different 
environments, will allow much more accurate measurements of the 
galaxy rotational velocity function in the local universe. 
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